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Using Monte Carlo techniques, the critical behaviour at edges and corners of the three- 
dimensional Ising model is studied. In particular, the critical exponent P2 of the local magnetization 
at edges formed by two intersecting free surfaces is estimated to be, as a function of the opening 
angle 9, 0.96 ± 0.02 for 6 = 135°, 1.28 ± 0.04 for 90°, and 2.30 ± 0.10 for 45°. The critical exponent 
/?3 of the corner magnetization of a cube is found to be 1.86 ± 0.06. The Monte Carlo estimates are 
compared to results of mean field theory, renormalization group calculations and high temperature 
series expansions. 
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I. INTRODUCTION 



Critical phenomena may be associated not only with the bulk of a three-dimensional crystal, but also with its 
surfaces, edges and corners. The critical properties of perfect, flat surfaces have been studied extensively [1-3]. 
However, the singularities of edge and corner quantities have attracted less attention. Some years ago, Cardy noted 
and calculated the dependence of the edge critical exponents on the opening angle 9 between the surfaces forming the 
y—i , edge (or wedge), using mean field theory and renormalization group theory of first order in e for 0{n) models [4]. Most 
of the later work dealt with Ising models in two dimensions [5] or polymers [6] . Rather rarely, the critical behaviour of 
edges in three-dimensional Ising models (describing magnetic crystals or lattice gas systems such as alloys) has been 
studied, applying renormalization group theory or high temperature series expansions [7-9]. Likewise, the singular 
properties at the corners of three-dimensional systems has been investigated even only in the framework of mean field 
O 1 tneor y I 5 ]- 

Edges and corners may be expected to play a dominant role, for instance, in nanostructured materials. Recently, 
Q\ , many studies have been performed to reproduce theoretically properties of such small clusters of atoms, including 
■ several Monte Carlo simulations [10,11]. In that context, it looks worthwhile to study systematically various types of 
^ \ edges and corners in somewhat simplified models, as well, as will be presented in this article. 

The outline of the article is as follows. In the next section, we shall consider three-dimensional Ising models with 
different opening angles, 9, at the edges. In particular, the critical exponent, /?2, of the edge magnetization is esti- 
mated, based on Monte Carlo (MC) data obtained from a cluster-flip algorithm. The values for are compared 
to previous results. In section 3, we shall present our Monte Carlo findings on Ising cubes, calculating especially the 
, corner magnetization and its critical exponent ft.3. Difficulties in extracting the bulk critical exponent from the total 
!LJ ' magnetization of the cubes are pointed out. A brief summary concludes the article. 



II. EDGES WITH DIFFERENT OPENING ANGLES 



We consider nearest-neighbour Ising models on simple cubic lattices with ferromagnetic interactions, J > 0. The 
Hamiltonian is given by 

(id) 

where the sum runs over all bonds («,j) between spins, Si — ±1, at neighbouring sites i = (x,y,z) and j. In the 
thermodynamic limit, the systems display a phase transition with the magnetization (S xyz ) vanishing above the critical 
point. From numerical analyses, the critical temperature, T c , is known very accurately, fcsT c = 4.511. .J [12,13]. 
To introduce edges in Ising magnets, we apply periodic boundary conditions along the, say, z-axis. The remaining 
four free surfaces of the crystal may be oriented in various ways leading to different opening angles 9 at the edges. In 
particular, for pairs of (100) and (010) surfaces, the opening angle is 9 — 90°, and the system contains four equivalent 
edges; case (a) in the following (see Figure 1). Obviously, the opening angle is invariant against rotating the crystal 
about the z-axis. However, for instance, the coordination numbers at the edges may change under a rotation, and 
therefore the local thermal edge properties may change. Keeping 9 = 90°, we rotated the crystal by 45° (yielding 
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(110) and (110) free surfaces), in order to check whether the critical behaviour at the edges is affected, case (6). In 
case (c), the surface orientations were chosen to be (100) and (110), with the intersections forming two pairs of edges 
with 9 — 45° and 9 — 135°. Note that edges with an opening angle of 180° correspond to the free surface, whose 
critical properties have been investigated rather extensively before [1,2,14,15]. 

In our present Monte Carlo study of Ising magnets with edges, we used the one-cluster-fiip algorithm [16]. Systems 
with L x M x N sites (or spins) were considered, with L x M being the number of sites in the planes perpendicular to 
the 2-axis, and N being the number of sites along that axis. Typically, L — M ranged from 20 to 60, with N going up 
to 640. The concrete values were chosen so that finite-size effects could be monitored and avoided, when attempting 
to elucidate critical properties of indefinitely extended edges, with the bulk and surface properties being those of 
the thermodynamic limit (see below). Usually, 4 x 10 4 clusters were generated per MC run, discarding the first 10 4 
clusters for equilibration. We averaged over an ensemble of at least five runs (using different random numbers) to 
obtain the final thermal averages. 

The exchange couplings are assumed to have always the same value, J > 0, i.e. in the bulk as well as at the surfaces 
and the edges. In the language of surface critical phenomena [1,2], one then encounters an ordinary transition, with 
the bulk, surfaces and edges ordering at the same (bulk) critical point, T c . Because one expects universal critical 
behaviour at the ordinary transition, that choice of couplings is merely for reasons of convenience. To cross over 
to different universality classes (corresponding to the special point or the extraordinary and surface transition), the 
surface couplings had to be increased significantly [1,2,14,15]. 
The crucial quantity of our analysis is the edge magnetization 
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where (x yoz) denotes a site at an edge of the lattice (in our notation, we follow Cardy [4], using the index '2' 
to indicate that an edge is formed by the intersection of two planes). In addition, we computed the analoguous 
magnetization for lines parallet to the z-axis, mi(x, y), as well as the susceptibility, \2, describing the response of the 
magnetization to a bulk field [4]. 

On approach to the bulk critical point, T c , m 2 is expected to vanish as m 2 oc t@ 2 , where t is the reduced temperature, 
t = \T C — T\/T c . To determine (32(9), we define, below and at the critical point, an effective exponent [14] as 

fo,eff(t) =dln{m 2 )/dln{t) (3) 

Because the MC data are recorded at discrete temperatures, U, (3 2 , e ff may be determined numerically by 

f3 2 ,eff(t)=Hm(z,t i )/m(z,t i+1 ))/ln{t i /t i+1 ) (4) 

with t = (ti+t i+1 )/2. 

02, e ff is well defined at all temperatures, not merely close to T c (for a related recent study of the corresponding 
critical exponent near criticality for ferromagnetic thin films, see Ref. [17]). Obviously, when t — ► 0, (3 2 , e ff becomes 
the asymptotic critical exponent (3 2 ■ In complete analogy, one may define the effective exponent of the line magne- 
tization, mi. For sufficiently large systems, m; approaches deep in the interior the bulk magnetization, nib, with the 
corresponding critical exponent (3b — 0.32.. [12,13]. At the surfaces, far from the edges, mi approaches the surface 
magnetization, mi, with the asymptotic critical exponent being (3\ « 0.80 [14]. 

To estimate reliably (3 2 from the simulational data, finite-size effects have to be monitored carefully. Strictly speaking, 
(3 2 is defined in the thermodynamic limit, L,M,N — ► oo (note that we take the absolute value to define the edge 
magnetization m 2 in (2), because we are simulating finite lattices). Of course, m 2 is affected by spin fluctuations in 
the bulk and at the surfaces. Accordingly, the MC system should be large enough to reproduce the thermodynamic 
values of nib and mi, sufficiently far away from the edge (these values are known very accurately [12-15]). In general, 
at a given temperature and geometry of the crystal, m 2 has to be stable against enlargening the system size, L x M 
and N . 

In Figure 2, the finite-size effect is illustrated for the geometry of case (c), with 9 = 45° and 135°, by plotting (3 2 , e ff 
versus the reduced temperature (error bars, resulting from ensemble averaging, were determined in the same way as 
before [14]). Finite-size dependences, i.e. deviations from the thermodynamic behaviour, are clearly signalled by a 
decrease of P 2 , e ff on approach to T c . The crucial quantity is evidently the length of the edge, N (in the range of 
temperatures depicted in the figure, the size of the planes perpendicular to the z-axis is sufficiently large, L = M = 40, 
to neglect further finite-size dependences, as we checked). By comparing our MC data for 9 = 45°, 90°, 135°, and 180° 
(free surface), an increasing edge length is needed to circumvent finite-size effects at fixed temperature and decreasing 
opening angle. 
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This trend may be understood by presuming that the critical behaviour at the edges is mainly governed by bulk 
critical fluctuations, similar to the situation at surfaces [14,15]. Now, by lowering 9, less pathes between neighbouring 
spins or sites connect the edges and the bulk. Accordingly, the edge length has to be enlarged to 'transmit' the bulk 
fluctuations fully to the edges at that reduced connectivity. 

It would be of interest to establish the finite-size scaling form of the edge magnetization m 2 (i, L = M,N) or, more 
generally, of the line magnetization mi(d,t,L = M,N), where d denotes the distance of the line to the edge. How- 
ever, such analyses are beyond the scope of this article.-It seems worth mentioning that our MC data suggest that 
the profile mi(d) approaches the surface magnetization, mi, or the bulk magnetization, m b , along, e.g., the shortest 
pathes, in an exponential form, i.e. 

m b (i) - mi(d) oc exp(-a b(1) d ). (5) 

From a preliminary analysis of our MC data, it seems conceivable that a b and a\ become identical on approach to T c , 
being related to the bulk correlation length. A similar behaviour holds for the magnetization profile as a function of 
the distance from the surface [14,15,18,19]. 

Figure 3 summarizes our simulational results for the effective exponent (32,eff(9,t) in the cases (a) and (c), i.e. with 
one pair of (100) surfaces. Only those MC data are shown which are not significantly affected by finite-size effects. 
We also included results of our previous MC study for the surface magnetization [14], corresponding to 8 = 180°. 
The slope of /3 2 , e // changes rather mildly at sufficiently small reduced temperatures t. Thereby meaningful esti- 
mates for the asymptotic exponent (32(0) seem to be feasible, with uncertainties rising when lowering the opening 
angle, because of a lack of MC data close to criticality due to the stronger finite-size effect. Our estimates arc 
fo{0 = 180°) = 0.80 ± 0.01 [14], /3 2 (135°) = 0.96 ± 0.02, /? 2 (90°) = 1.28 ± 0.04, and /3 2 (45°) = 2.30 ± 0.10, with error 
bars referring to 'reasonable' extrapolations of the effective exponents to the critical point. 

These MC estimates are to be compared with results obtained from mean field theory (MF) [4] , renormalization group 
calculations to first order in e (RNG) [4], and high temperature series expansions (HTS) [7], see Figure 3 and Table 
1. The values given by mean field theory are presumably systematically too high, according to the predictions of the 
renormalization group. These predictions, in turn, seem to be systematically too large, as suggested both by high 
temperature series expansions and the MC simulations. The last two methods yield results which are in close agree- 
ment with each other. A similar tendency in the values of /3 2 (#) obtained from the different analytical and numerical 
approaches (MF, RNG, HTS, and MC) holds for the n-vector model with n = 0, describing polymers [4,7,20]. 
At the fixed opening angle 8 = 90°, we studied the effect of a rotation of the crystal on the critical edge properties, 
case (6). In particular, the crystal was rotated by 45° about the z-axis. The wedge is then formed by (110) and (110) 
surfaces. By the rotation, the coordination numbers at the surfaces and edges are reduced by one. Consequently, the 
surface and edge magnetizations are lowered. Following the above considerations, the finite-size effect is expected 
to be enhanced compared to case (a), i.e. after rotation, larger MC systems are needed to approximate closely the 
thermodynamic behaviour. The MC data confirm this expectation. 

Results on the effective exponent /3 2 , e // are displayed in Figure 4, comparing the cases (a) and (6). Again, we depict 
only MC data which are not significantly affected by finite-size effects. It seems well conceivable to extrapolate both 
cases to the same asymptotic critical exponent, /3 2 = 1.28 ± 0.04, consistent with the invariance of the value of the 
boundary critical exponent against rotation of the crystal. After rotation, we did not approach the critical point as 
closely as before, because otherwise extremely large systems had to be dealt with. In any event, as seen in Figure 4, 
corrections to scaling are much weaker in case (6), for small reduced temperatures, resulting in the much lower slope 
of /3 2>e //, compared to Ising models with (100) and (010) surfaces. 

Note that the invariance of boundary critical exponents against rotation of the crystal has been discussed and partly 
even proven for surfaces and edges of two-dimensional Ising models [5,19,21]. 

Our MC data for the exponent of the edge susceptibility, xi-> are l ess accurate than those for the edge magnetization, 
thereby allowing to estimate the corresponding critical exponent 7 2 only with rather large uncertainties. At any rate, 
7 2 is found to depend strongly on the opening angle 9, and the results of the RNG calculations and HTS method [4,7] 
seem to provide fairly good estimates for the true values. 



III. CUBES AND CORNERS 

We consider Ising cubes of L 3 spins on simple cubic lattices with free surfaces. In the thermodynamic limit, 
L — > oo, the corner magnetization is given by 
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Sx a y a z a ) (6) 

summing and averaging over the eight corner spins at site (xoj/o^o)- I n complete analogy to (3), one may define an 
effective exponent /33. e //, approaching the asymptotic corner critical exponent /?3, as t — > 0. 

In a finite system, m 3 is, strictly speaking, zero, because the Hamiltonian is invariant against reversing the spins. On 
the other hand, for instance, the absolute corner magnetization 
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seems to show no critical behaviour. Indeed, 7773. a is quite distinct from 7713, involving a fixed number of spins 
independent of the system size, in contrast to the bulk, surface or edge magnetizations calculated with or without 
absolute values. 

To compute m 3 , we use the standard one-spin-flip MC algorithm (Metropolis algorithm). Taking the fully ordered 
ground state with, say, S xyz — 1 at all sites as starting configuration of the simulation, one encounters, after initial 
relaxation, at T < T c , a metastable state, in which the total magnetization is positive. The system remains in 
that metastable state (being separated by an energy barrier from the 'mirror' metastable state with negative total 
magnetization), possibly, for a long time, the time depending on the size L and the temperature, see, e.g., Ref. [22]. 
The corresponding local magnetizations, especially the corner magnetization 777,3, do not vanish then, and their values, 
taking usual finite-size dependences into account, are supposedly very close to those in the thermodynamic limit. 
In fact, by monitoring the total magnetization during the MC runs, we identified the temperatures at which we can 
reliably compute 7773. In Figure 5, resulting findings for the effective exponent, /?3 )e //, are shown. Most of the data 
are for rather small cubes (mimicking magnetic nanoparticles) , with L = 20. 3000 Monte Carlo steps per site (MCS) 
were used to 'equilibrate' (to the metastable state) the systems; additional 7000 MCS were performed to obtain 
thermal averages. To improve the statistics, we sampled over an ensemble of a few 10 2 independent MC runs at each 
temperature. No finite-size effects were observed, when we increased L to 40 at the smallest reduced temperatures, 
see Figure 5. To generate data of the high accuracy required to determine reliably the effective exponent without 
finite-size corrections closer to T Cl one had to simulate eventually much larger system sizes. In that way, the computing 
time had to be increased appreciably then. We refrained from doing so, because the present data already allow for a 
decent estimate of the asymptotic critical exponent, (3% = 1.86 ± 0.06. That value is significantly lower than the one 
predicted by mean field theory, (3^ — 2 [5] . 

In finite-size scaling theories [23,24], one usually considers the total magnetization of the system, summing and 
averaging over all spins. This approach may work in the simplest form by using periodic boundary conditions, implying 
translational invariance. Of course, in the case of free surfaces much care is needed, because the contributions of bulk, 
surface, edge, and corner spins to the total magnetization arc of different type. Actually, the local magnetization 
differs, in general, from site to site. 

In Figure 6, the effective exponent, /?T,e// of the total magnetization tot, taking, in analogy to (2), the absolute 
values, is shown again for a relatively small cube with 20 3 spins (including, for comparison, data for L — 40). (3r,eff 
is appreciably higher than the bulk critical exponent (3b = 0.32.. already quite far away from criticality, t < 0.2. 
The 'overshooting phenomenon' is caused by the spins near the surfaces, edges and corners, reducing the total 
magnetization and yielding large effective and critical exponents, as discussed in this article. Certainly, on approach 
to T c the effective exponent goes to zero, because the absolute total magnetization tut does not vanish at the critical 
point in a finite system. 

The overshooting persists at larger system sizes (becoming, however, weaker at fixed t; see Figure 6 for L =40), and 
it may be difficult to extract the correct bulk critical exponent reliably In an alternate approach, one may replace 
the reduced temperature t by t m = \T — T m \/T m , where T m denotes the turning point of the total magnetization for 
a given system size, L. From (3 e ff = efln (m^/dln (t m ), one might estimate (3b (for L — > 00 one has t m — > t, and 
thence /3 e // — ► (3b as t — ► 0) in an easy way, see Figure 6. For larger system sizes, (3r,eff and (3 e ff tend to approach 
each other, as one observes, e.g., for systems with L — 40. In any event, detailed considerations would be desirable 
to substantiate and quantify the finite-size effects for Ising cubes (for related problems, demonstrating the sensitivity 
of finite-size scaling to different boundary conditions, see, e.g., Ref. [23,24]). 

IV. SUMMARY 

Using Monte Carlo techniques, we studied critical properties of edges and corners in three-dimensional Ising models. 
In particular, by computing effective exponents for the edge and corner magnetization, the corresponding asymptotic 
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critical exponents could be estimated quite reliably. 

The critical exponents at the edges are found to depend sensitively on the opening angle formed by the surfaces 
intersecting at the edge. The numerical values for these exponents, based on the simulational data, refined predictions 
of renormalization group calculations of first order in e, and confirmed results of high temperature series expansions. 
The exponents seem to be invariant against rotating the crystal, but keeping the opening angle, as we demonstrated 
in a specific case. However, the corrections to scaling may be largely affected by such a rotation. 
The critical exponent at the corner had been calculated before only by using mean field theory. The estimate based 
on the Monte Carlo data shows that the true value is significantly lower than the one predicted by mean field theory. 
It would be quite interesting to apply other numerical and analytical techniques to determine that exponent. 
By simulating corner properties, we considered Ising models on a cube with free surfaces (mimicking magnetic nanopar- 
ticles). Such systems show very strong finite size effects, when calculating quantities averaging over the entire lattice. 
Care is needed in extracting bulk critical properties, because the singularities at the surfaces, edges, and corners are 
distinct from the bulk singularities, and play an important role even for fairly large system sizes. 
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TABLE I. Predictions for edge critical exponent P2, using various methods discussed in the text. 





45° 


90° 


135° 


180° 


MF [4] 


2.50 


1.50 


1.17 


1.00 


RNG [4] 


2.48 


1.39 


1.02 


0.84 


HTS [7] 


2.30 


1.31 


0.98 


0.81 


MC 


2.30 ±0.10 


1.28 ±0.04 


0.96 ±0.02 


0.80 ± 0.01 [14] 



Figure Captions 

Fig. 1: Geometry of an Ising model with (100) and (010) surfaces (shadowed), i.e. edges with opening angle 9 = 90°. 

Fig. 2: Effective exponent (32, e ff of the edge magnetization versus reduced temperature, t, at edges with opening angles 
9 = (a) 45° and (b) 135°, simulating different sizes L x M x N, see inset. 

Fig. 3: fi2,eff versus t at edges with opening angles 9 — 180° (triangles) [14], 135° (diamonds), 90° (squares), and 45° 
(circles). MC system sizes were chosen to circumvent finite-size effects, see text. The predictions of RNG 
calculations (dashed line) [4] and HTS expansions (dot-dashed line) [7] for /3 2 are shown as well. 

Fig. 4: fl2,eff versus t at edges with 9 = 90°, with pairs of (100) and (010) surfaces (open circles) and pairs of (110) 
and (110) surfaces (full circles), for MC system sizes up to 60 x 60 x 640 to avoid finite-size dependences. 

Fig. 5: Effective exponent fi^^ff of the corner magnetization as a function of t for Ising cubes with 20 3 spins (open 
circles) and, for comparison at small t, with 40 3 spins (full triangles). The broken line is a guide to the eye. The 
solid line refers to the value of /?3 as calculated from mean field theory [5]. 

Fig. 6: Effective exponent of the total magnetization f3r, e ff (Peff', full symbols) versus reduced temperature t (t m ) of an 
Ising cube with 20 3 (squares) and 40 3 (circles) spins, using different reference temperatures for the (finite-size) 
critical point T c (T m ), see text. 
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